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1.  RESEARCH  OBJECTIVE 


The  overall  objective  of  the  proposed  research  is  to  develop  and  verify 
mathematical  models  of  deformation  and  delamination  of  elastic  and 
viscoelastic  fibrous  composites  with  distributed  damage.  Emphasis  of  the 
theoretical  and  experimental  studies  is  on  graphite/epoxy  composites  under 
complex  loadings  which  produce  matrix  damage  and  delamination. 


2.  ACCOMPLISHMENTS 


2.1  Overview  . -.-.L 

Methods  of  deformation  and  fracture  characterization  are  simplified  when 
strain  energy-like  potentials  based  on  mechanical  work  can  be  used.  Research 
during  the  grant  period  has  been  primarily  concerned  with  development  of  the 
theoretical  basis  for  the  approach,  and  with  experimental  studies  which 
demonstrate  its  validity  for  rubber-toughened  and  untoughened  granhi te/epoxy 
composites  with  growing  or  constant  matrix  damage  and  delamination. 

In  the  following  sections  the  research  work  and  primary  findings  are 
summarized  by  abstracting  the  papers  and  dissertations  which  detail  the 
studies.  In  Section  2.2  the  general  theory  is  covered;  it  allows  for  evolving 
or  constant  damage.  Viscoelastic  effects  in  the  matrix  with  constant  damage 
are  considered  in  Section  2.3  in  a  study  of  tubes  under  axial  and  torsional 
loading.  Sect"'""  2.4  discusses  predictions  of  matrix  damage  due  to  thermal 
stresses  and  an  .-rimental  program  that  supports  the  theory.  An  extensive 
experimental  study  of  bars  under  axial  and  torsional  loading  is  described  in 
Section  2.5;  path- independence  of  work  during  periods  of  growing  damage  is 
shown  for  composites  with  toughened  and  untoughened  resins.  Sections  2.6  and 
2.7  cover  predictions  and  experimental  verification  for  in-plane  loading  of 
laminates;  the  method  of  analysis,  which  is  based  on  a  work-potential  and  its 


2 


minimization,  is  compared  in  the  former  study  to  a  method  of  analysis  based  on 
plasticity  theory.  The  work-potential  provides  the  basis  for  the  J-integral 
in  delamination  fracture  analysis.  Investigations  using  this  type  of  analysis 
and  experimental  studies  of  the  double  cantilever  beam  with  multidirectional 
plys  is  summarized  in  Section  2.8.  Two  related  studies  of  beam-width  effects 
and  shear  deformation  effects  are  described  in  Sections  2.9  and  2.10, 
respecti vely. 

The  Appendix  contains  the  two  papers  which  appeared  after  the  annual 
report  was  submitted  in  August  1988. 

2.2  A  Theory  of  Mechanical  Behavior  of  Elastic  Media  with  Growing  Damage 
and  Other  Changes  in  Structure  [1]* 

Strain  energy-like  potentials  are  used  to  model  the  mechanical  behavior 
of  linear  and  nonlinear  elastic  media  with  changing  structure,  such  as  micro- 
and  macrocrack  growth  in  monolithic  and  composite  materials.  Theory  and 
experiment  show  that  the  applied  work  for  processes  in  whicn  changes  in 

structure  occur  is  in  certain  cases  independent  of  some  of  the  deformation 
history.  Consequences  of  this  limited  path-independence  are  investigated,  and 
various  relationships  for  stable  mechanical  response  are  derived.  For 

example,  it  is  shown  that  work  is  a  minimum  during  stable  changes  in 

structure,  which  should  be  useful  for  developing  approximate  solutions  by 
variational  methods.  Some  final  remarks  indicate  how  the  theory  may  be 

extended  to  include  thermal,  viscoelastic,  and  fatigue  effects. 

2.3  A  Method  For  Studying  Composites  with  Changing  Damage  by  Correcting 
For  the  Effects  of  Matrix  Viscoelasticity  | 2) 

A  technique  is  described  for  modifying  stress-strain  data  on  fibrous 
composites  so  that  effects  of  changing  damage  may  be  observed  without  the 

♦The  number  in  brackets  refers  to  the  publication  which  is  abstracted  here; 
the  publication  list  starts  on  page  6. 
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complicating  effects  of  matrix  viscoelasticity.  The  method,  which  is  based  on 
micromechanical  considerations,  reduces  the  behavior  to  that  of  an  equivalent 
elastic  composite  with  damage.  The  fibers  are  assumed  to  be  continuous  and 
linearly  elastic.  The  theoretical  basis  is  developed  and  then  the  method  is 
illustrated  using  results  from  cyclic  axial-torsional  loading  of  tubular 
specimens  of  graphite/epoxy  laminates. 

2 . 4  Thermal ly- Induced  Fracture  in  Composites  [3] 

This  work  is  an  experimental  and  analytical  investigation  of  thermally- 
induced  cracking  in  cross-ply  graphite/epoxy  composites.  It  is  shown 
experimentally  that  both  rates  and  amplitudes  of  the  thermal  excursions  alf'ct 
the  extent  and  the  form  of  damage.  The  analytical  study  shows  that  the  early 
stages  of  sufficiently  slow  thermal  excursions  result  in  crack  patterns  that 
are  analogous  to  mechanical  loading  effects,  and  can  be  assessed  by  an 
approximate,  two-dimensional  micro-cracking  model.  However,  three-dimensional 
aspects  of  the  spatial ly  non-uniform  stress  field  may  have  to  be  included  to 
model  crack  formation  under  subsequent  temperature  excursions  or  rapid  thermal 
fluctuations.  In  the  latter  cases  oblique  and  curved  cracks  develop  and  the 
laminate  is  susceptible  to  internal  and  free-edge  delaminations. 

2 . 5  Deformation  and  Delamination  of  Inelastic  Laminates  Under  Tensile 

and  Torsional  Loading  [4,5] 

This  study  is  a  fundamental  examination  of  the  theoretical  hypothesis 
that  mechanical  work  is  a  (multivalued)  potential  function  which  characterizes 
the  deformation  and  fracture  behavior  of  inelastic  materials  during  damage 
growth  processes.  Experimental  data  from  fiber-reinforced  plastic  laminates 
(with  and  without  rubber  particle  toughening)  subjected  to  axial  and  torsional 
deformation  are  analyzed  for  the  existence  of  the  work  potential.  The  work 
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potential  is  first  employed  to  analyze  data  from  proportional  deformation 
tests  and  make  predictions  of  load  response.  Good  agreement  is  obtained 
between  theory  and  experiment.  Data  from  nonproportional  deformation  tests 
are  then  evaluated  for  displacement-path  independence  of  work  and  load 
response.  Domains  of  path  independence  are  found  from  specimens  strained  well 
into  the  range  of  nonlinear  inelastic  behavior.  Thus,  the  results  of  these 
experimental  studies  support  the  existence  of  a  work  potential.  Although 
viscoelastic  effects  are  present,  they  are  minimized  by  using  isochronal  data 
in  the  characteri zati on  of  mechanical  work.  Finally,  work  potential  theory 
and  experimental  results  are  used  to  determine  critical  energy  release  rates 
for  mixed  mode  delamination  of  laminates  subjected  to  axial  and  torsional 
defoi mat  ions . 

2 . 6  A  Method  for  Mechanical  State  Characterization  of  Inelastic  Composite 

Laminates  with  Damage  [ 6 ] 

The  method  using  a  work  potential,  and  its  minimization,  is  described  for 
the  characterization  of  mechanical  behavior  of  inelastic  composites  with 
damage,  but  without  significant  time-dependent  behavior.  It  is  based  on  the 
theoretically  and  experimentally  motivated  assumption  of  path-indepenuence  of 
mechanical  work  over  limited  ranges  of  stress  or  strain  states.  This  method 
and,  for  comparison,  an  approach  employing  plasticity  theory  are  illustrated 
with  the  special  case  of  a  unidirectional-fiber  laminate  or  ply.  Use  of  the 
work-potential  method  for  a  multidirectional-fiber  laminate  is  discussed  in 
the  concluding  remarks. 

2.7  Mechanical  Characterization  and  Analysis  of  Inelastic  Composite 

Laminates  with  Growing  Damage  [7) 

A  method  of  laminate  characterization  and  analysis  is  described  in  which 
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growing  damage  and  other  inelastic  phenomena  are  treated  using  the  same 
mathematical  formalism,  thus  simplifying  the  description  of  mechanical 
response.  It  is  based  on  the  observation  that  the  applied  work  is  not 
sensitive  to  many  details  of  the  deformation  history.  Following  a  brief 
discussion  of  the  thermodynamically-based  theory,  a  special  version  is  used 
along  with  experimental  data  on  graphite/epoxy  composite  to  obtain  an  explicit 
mathematical  characterization  of  a  unidirectional  ply.  Predictions  of 
mechanical  response  are  then  compared  to  experimental  results  for  a  variety  of 
layups,  one  of  which  delaminates  from  the  edges.  Good  agreement  between 
theory  and  experiment  is  shown. 

2 . 8  Determination  of  the  Mode  I  Delamination  Fracture  Toughness  of 

Multidirectional  Composite  Laminates  [8] 

The  objective  of  this  study  is  to  develop  and  verify  a  J-integral  method 
for  characterizing  mode  I  delamination  fracture  of  composite  laminates  with 
distributed  matrix  damage.  Attention  focuses  on  the  special  problems 
associated  with  delamination  of  composites  with  multidirectional  (as  opposed 
to  unidirectional)  layups.  Nonlinear  beam  theory  is  used  to  analyze  the 
double  cantilever  beam  specimen  to  derive  an  approximate  expression  for  the  J- 
integral.  A  related  test  method  is  proposed.  An  experimental  program  and 
results  are  described  which  explore  the  utility  of  the  method  and  the 
variables  affecting  delamination  of  multidirectional  composites. 

2 . 9  Effect  of  Finite  Width  on  Deflection  and  Energy  Releases  Rate 

of  an  Orthotropic  Double  Cantilever  Specimen  [91 


The  problem  of  an  orthotropic  cantilevered  plate  subjected  to  a  uniformly 
distributed  end  load  is  solved  by  the  Rayleigh-Ritz  energy  method.  The  result 
is  applied  to  laminated  composite,  double  cantilevered  specimens  to  estimate 
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the  effect  of  crack  tip  constraint  on  the  transverse  curvature,  deflection  and 
energy  release  rate.  The  solution  is  also  utilized  to  determine  finite  width 
correction  factors  for  fracture  energy  characterization  tests  in  which  neither 
plane  stress  nor  plane  strain  conditions  apply. 

2.10  A  Technique  for  Predicting  Mode  I  Energy  Release  Races  Using  A 
First  Order  Shear  Deformable  Plate  Theory  [10] 

Utilizing  a  first  order  shear  deformable  plate  theory,  a  technique  is 
described  for  predicting  the  distribution  of  the  energy  release  rate  along  a 
curved  or  straight  mode  I  planar  crack  in  the  plane  of  a  plate  (such  as  a 
delamination  crack).  Accuracy  of  the  technique  is  assessed  by  comparing  the 
distributions  of  energy  release  rate  with  those  predicted  by  two  and  three 
dimensional  finite  element  analyses  of  double  cantilever  beam  specimens  with 
straight  crack  fronts. 
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ABSTRACT 

A  method  using  a  work  potential  is  described  for  the  characterization  of 
mechanical  behavior  of  inelastic  composites  with  damage,  but  without 
significant  time-dependent  behavior.  It  is  based  on  the  theoretically  and 
experimentally  motivated  assumption  of  path-independence  of  mechanical  work 
over  limited  ranges  of  stress  or  strain  states.  This  method  and,  for 
comparison,  an  approach  employing  plasticity  theory  are  illustrated  with 
the  special  case  of  a  unidirectional-fiber  laminate  or  ply.  Use  of  the 
work-potential  method  for  a  multidirectional-fiber  laminate  is  discussed  in 
the  concluding  remarks. 


KEYWORDS 


Composites,  Laminates,  Damage,  Inelasticity,  Plasticity 
1.  INTRODUCTION 


Considerable  progress  has  been  made  in  recent  years  on  the  development  of 
high  strength-to-weight ,  tough  structural  composites.  This  behavior  is 
achieved  in-part  by  laminating  individual  plies  of  unidirectional, 
continuous  fiber-reinforced  plastic  or  metal.  The  laminates  are  resistant 
to  crack  growth  through  the  thickness  if  two  or  more  fiber  orientations  are 
used.  Delamination  and  cracking  within  plies  is  reduced  by  using  ductile 
matrices.  For  organic  polymer  matrices,  the  ductility  is  obtained  by 
adding  toughening  agents,  such  as  rubber  particles,  to  normally  brittle 
crosslinked  resins,  or  by  using  resins  with  little  or  no  crosslinking 
(Johnston,  1987).  These  improvements  in  material  performance  place 
increased  demands  on  the  structural  designer  and  those  concerned  with  the 
micromechanics  of  composites  if  inelasticity  is  due  to  both  plastic 
deformation  and  damage  or  if  it  has  to  be  considered  under  a  wider  range  of 
conditions  than  for  the  brittle  matrix  systems. 


Traditiona 

plasticity 

composites 


ly,  matrix  ductility  has  been  treated  using  incremental 
theory  (Christensen,  1979)  while  micro-  and  macrocracking  of 
have  been  analyzed  using  linear  elasticity  theory  (Wang  and 
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Haritos,  1987).  In  this  paper  we  discuss  an  approach  to  characterizing 
inelastic  composite  material  behavior  which  is  based  on  total  rather 
than  incremental  strains.  Also,  the  approach  uses  the  same  mathematical 
formalism  for  inelasticity  due  to  plastic  deformation  as  due  to  cracking  on 
various  scales  and  other  damage  mechanisms;  the  term  inelasticity,  as  used 
here,  refers  to  any  stable  behavior  in  which  stress  or  load  is  not  always  a 
single-valued  function  of  strain  or  displacement.  It  is  believed  that  this 
unified  approach  simplifies  the  problem  of  understanding  and  predicting 
mechanical  behavior  of  composites  with  damage.  Fatigue  and  time-dependent 
behavior  and  thermal  effects  are  not  treated  here,  although  approaches  have 
been  proposed  in  the  papers  which  motivated  the  present  study  (Schapery, 
1987a,  1988). 

Schapery  (1987a)  has  shown  theoretically  that  the  stresses  and  mechanical 
work  of  deformation  are  often  independent  of  many  details  of  the 
deformation  history  when  the  inelasticity  is  due  to  micro-  and 
macrocracking.  However,  cracking  is  not  the  only  mechanism  that  produces 
this  behavior.  Indeed,  it  has  been  observed  for  a  rubber-toughened, 
graphite/epoxy  composite  in  which  there  are  probably  significant  effects  of 
shear  banding  in  the  matrix  (possibly  initiated  or  enhanced  by  cavitation 
of  rubber  particles)  (Yee,  1987).  This  limited  path-independence  was  used 
by  Schapery  (1988)  to  develop  a  constitutive  theory  that  treats  different 
inelastic  mechanisms  within  the  same  mathematical  framework.  Also,  as 
shown  by  Schapery  (1987a),  fracture  analysis  is  simplified  when  this  theory 
is  valid  because  of  the  applicability  of  certain  equations  for  relating 
changes  in  local  and  global  energies. 

Figure  1  illustrates  one  type  of  path-independence  we  have  found  for  the 
rubber-toughened  composite.  Rectangular  composite  bars  with  an  angle-ply 
layup  (alternating  fiber  angle,  e  =  t  35°,  with  respect  to  the  axial 
direction)  were  subjected  to  various  axial  and  torsiona1  deformations 


Fig.  1.  Shear  stress-strain  curves  for  proportional  and  nonpro¬ 
portional  straining  of  an  angle-ply  laminate;  Hexcel  T2C 
145/F155  graphite/epoxy  [ +- 35 0  ] gs ;  0.15"  thick  X  0.5" 
wide  X  8.75"  long.  From  Lamborn  and  Schapery  (1988). 
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through  controlled  movement  of  the  end-grips.  The  different  deformation 
paths  are  identified  in  Fig.  1  by  number;  for  example,  the  bottom  line  type 
is  used  for  axial  history  1  during  the  first  loading  period  and  axial 
history  6  during  unloading,  while  the  corresponding  torsional  histories  are 
2  and  5.  The  "nominal"  shear  stress  and  shea'"  and  axial  strains  are 
quantities  which  are  proportional  to  the  torque,  twist,  and  axial 
displacement,  respectively;  the  proportionality  coefficients  depend  only  on 
the  specimen  dimensions,  and  are  introduced  to  minimize  the  effect  of 
specimen-to-specimen  size  differences. 

At  the  end  of  the  first  loading  period,  the  five  different  strain  paths 
result  in  practically  the  same  stress  (Fig.  1)  and  total  work.  The  same 
behavior  holds  for  the  unloading  and  reloading.  In  contrast,  unreinforced 
aluminum  bars  exhibit  significant  path-dependence  (Lamborn  and  Schapery, 
1988);  we  do  not  know  if  f i ber-reinforced  aluminum  would  exhibit  less  path- 
dependence. 

Unloading  and  reloading  behavior  of  the  graphite/epoxy  material  under  pure 
axial  or  torsional  straining  is  similar  to  that  shown  in  Fig.  1;  there  is 
significant  hysteresis  and  the  average  slope  of  the  loop  decreases  with 
increasing  strain  at  the  unloading  point.  The  stress  during  loading  does 
not  usually  exhibit  a  maximum  point  prior  to  fracture,  in  contrast  to  that 
in  Fig.  1.  We  are  now  investigating  the  damage  state  as  a  function  of 
deformation  history  using  similar  specimens;  significant  edge  delaminations 
have  been  found  at  the  highest  stresses  for  deformation  histories  like 
those  in  Fig.  1. 

The  primary  effects  of  deformation  history  on  the  composite  appear  to  be 
associated  with  the  sign  of  (nominal)  strain  rate  and  the  strain  magnitude 
when  the  sign  last  changed.  Although  a  more  precise  definition  of  limited 
path-independence  was  given  by  Schapery  (1988)  here  we  shall  just  refer  to 
differences  between  loading,  unloading,  and  reloading  curves,  and  suppose 
that  for  each  case  there  is  no  effect  of  path  (which  is  approximately  true 
for  the  data  in  Fig.  1). 

The  local  stresses  and  strains  (as  opposed  to  the  "nominal"  quantities  in 
Fig.  1)  are  distributed  very  nonuniformly  throughout  the  specimens  used  in 
these  axial-torsional  tests,  and  thus  the  results  cannot  be  used  directly 
in  a  basic  material  characterization  of  the  composite.  However,  it  is 
unlikely  that  the  specimens'  overall  behavior  would  exhibit  limited  path- 
independence  if  the  ply-level  constitutive  equations  did  not  reflect  this 
type  of  behavior. 

The  discussion  in  Sections  2-4  is  concerned  primarily  with  the 
characterization  of  the  behavior  of  a  undirectional -fiber  laminate 
consisting  of  one  or  more  plies  under  the  assumption  of  this  limited  path- 
independence.  Special  versions  of  the  theory  (Schapery,  1988)  are  used  here 
to  illustrate  it  for  composites.  Specifically,  Section  2  considers 
nonlinear  loading  and  unloading  behavior,  and  expresses  the  inelasticity  in 
terms  of  one  parameter  S  which  represents  the  effect  of  microstructural 
changes  on  the  overall  stress-strain  behavior;  such  S-parameters  provide 
the  inelasticity  and,  in  the  context  of  some  thermodynamic  formulations, 
are  called  internal  state  variables.  Section  3  contrasts  the  theory  with  a 
plasticity  model  based  on  the  normality  rule,  and  uses  the  characterization 
in  Section  2  as  an  example.  In  Section  4  another  illustration  is  given  by 
using  a  linear  approximation  for  unloading  behavior.  Concluding  remarks  in 
Section  5  discuss  in-part  the  use  of  unidirectional  ply  characterization  in 
laminates  with  ply-level  and  larger  scales  of  damage. 
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2.  A  CONSTITUTIVE  EQUATION  WITH  NONLINEAR  UNLOADING  BEHAVIOR 

Figure  2  shows  a  unidirectional  laminate  or  ply  and  the  coordinate 
notation,  in  which  the  axis  is  parallel  to  the  fibers;  the  x^  axis  is 

normal  to  the  ply  plane.  The  stresses  o.  and  strains  (i  =  1,2, ...6) 

are  mechanical  variables  referred  to  the  principal  material  coordinates 
x^ .  In  most  of  the  discussion  it  will  be  convenient  to  use  this  single 
index  notation.  As  is  customary,  i  =  4,5,6  are  used  for  the  shearing 
variables;  the  relationship  between  single  and  double  indexed  variables  for 
plane  stress  is 


«— H 

D 

II 

O 

°22  ~  °2* 

°12'  °6 

(1) 

£11  "  £1* 

e22  =  e2  * 

2£12  =  £6 

A  constitutive  equation 

will  be 

proposed  which 

accounts  for  nonlinear 

loading  and  unloading  behavior  and  which  is  consistent  with  the  path- 
independence  of  work  discussed  in  the  Introduction  as  well  as  the  nonlinear 
behavior  reported  by  Lou  and  Schapery  (1971)  and  Sun  and  Chen  (1987);  the 
reader  is  referred  to  these  two  papers  for  the  experimental  data,  as  space 
does  not  permit  its  reproduction  here.  Specifically,  a  strain  °nergy 
density  w  =  w(e.,  S)  is  assumed  to  exist,  where  the  microstructure  state  is 
defined  by  S;  only  one  structure  parameter  S  will  be  used  here,  although 
more  could  be  introduced,  if  necessary.  By  definition  of  w. 


oi  =  3w/ 3  c i  (2) 

In  both  aforementioned  references  strains  are  expressed  in  terms  of 
stresses,  and  thus  it  is  helpful  to  eliminate  w  in  favor  of  a  so-called 
dual  strain  energy  density  wq  =  w^o^  S) , 

w0  E  w  -  0^  (3) 

(Throughout  this  paper  the  summation  convention  is  employed,  in  which  a 
repeated  index  implies  summation  over  its  range.)  By  using  (2)  and 
introducing  differential  changes  in  (3),  it  follows  in  the  usual  way  that 


e,=  -  3W  / 3 o . 
1  0  1 


(4) 


A  form  of  w  discussed  by  Schapery  (1988,  Eq.  (A24))  is  proposed  now  for 
characterizing  ply  behavior. 
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Fig.  2.  Unidirectional  composite  and  coordinates. 


5 


w  =  w  +  P(c  .  S) 
o  cO  '  0 


where  w  =  (o . ) ,  o  =  a  (a.),  and  P  are  presently  arbitrary  func¬ 
tions.  °?ne  nfe°chanicar  worR  during  processes  in  which  S  changes  can  be 
shown  to  be  independent  of  path  if  and  only  if  S  =  S(o  ).  ;  proof  of  this 
statement  may  be  made  by  the  same  method  as  used  in  a  study  of  w  (Schapery, 
1988,  Appendix  A).  The  function  S(o  )  can  be  absorbed  in  the  functional 
dependence  of  P  on  S,  and  thus  we  may  use  S  =  a  whenever  S  changes  without 
any  actual  limitation  in  the  model.  Whether  S° varies  or  is  constant,  the 
strains  are  obtained  from  (4)  and  (5), 


e  aP  3°o 

£  i  C  i  3  a  3  a  ■ 

0  1 


where,  by  definition, 


£i  5  ‘  3woo/3oi  (7) 

p 

The  e-  are  defined  through  derivatives  of  a  fully  path-independent 
potential,  w  ,  and  thus  it  is  appropriate  to  call  them  "elastic  strains". 
All  stress-hS^tory  effects  are  in  the  second  term  in  (6),  which  gives  the 
"inelastic  strains". 

In  order  to  obtain  a  constitutive  equation  that  agrees  with  Sun  and  Chen's 
experimental  data  we  select  for  oQ  the  quadratic  form, 

°o  =  (aij  °i  (8) 

where  the  a ^ j  are  constants;  as  the  antisymmetric  components  of  a^  have  no 
effect  on  o_  ,  we  may  suppose  a^-  =  a^.  Now, 


3 °o/3 °i  =  aij  °j  /oo 


and  thus  from  (6), 


e  3P  , 

£  .  =  £  .  -  -  a  •  •  C  .  /  O 

1  1  3o0  1J  J  0 


During  structure-change  processes  S  =  o  ,  as  noted  previously,  and 
therefore  the  coefficient  3P/3oQ  depends  on°only  oQ.  For  such  processes 
we  may  thus  write 


where 


e .  +  e  a .  •  o  •  /o 
i  o  ij  j  '  o 


e  =  e  ( o  )  r  -(3P/3o„)  evaluated  at  S  =  o„ 
0  o  0  '  0  0 


In  the  terminology  of  plasticity  theory,  (11)  is  for  "loading"  processes. 

Without  fiber  fracture,  the  strain  in  the  fiber  direction  is  essentially 
independent  of  stress  history  in  most  structural  composites;  thus,  as 
assumed  by  Sun  and  Chen,  a,  .  =  a  .,  =  0.  We  suppose  further  that  the 
composite  is  orthotropic,  regard lesrW  stress-history,  where  the  axes  xi 
are  the  principal  material  axes;  this  condition  implies  the  only  a^  which 

do  not  vanish  are  &22'  a23’  ^33’  a44,a55»  a66*  W^1  T  3S  a. 3^ ( =  ^2 3 ^ 

are  really  only  five  independent  constants  because  o  may  be  normalized 
with  respect  to  a  constant  without  limiting  the  generality  of  (5);  this 
normalization  will  be  done  by  simply  letting  &22  =  If  a11  stresses 
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vanish  except  for  o~,  (8)  reduces  to  a  =  |aJ;  thus  o  becomes  the 
applied  stress  for  tire  case  of  uniaxial  Pensile  loading  normal  to  the 
fibers . 

For  plane  stress,  °3  =  °4  =  °5  =  so  that  (8)  reduces  to 


co  =  (o22+  a66  °62)H 


From  (11), 


e 

el  =  £1 


£2  =  c2  +  co°2/oo  (15) 

£6  =  £6  +  a66  eo°6/oo  (16) 

For  uniaxial  tension  normal  to  the  fibers,  o  =  o~,  as  noted  previously. 
Equation  (15)  then  shows  that  e  reduces  to  °the  inelastic  component  of 
e9.  For  general  stress  states  e°  is  at  most  a  function  of  o  ,  according 
to^( 12) .  0  0 

By  introducing  some  additional  specializations,  including  the  assumption 
that  the  e ^  are  linear  in  the  o. ,  we  will  finally  arrive  at  Sun  and  Chen's 
findings  for  uniaxial  loading  of  unidirectional,  rectangular  specimens. 
Namely,  for  loading  in  the  x  direction  (cf.  Fig.  2), 

2  .2 

0^  =  cos  6  g^,  ~  s^n  «  °x»  =  -sine  COSeo^  (17) 

where  a  is  the  applied  force/area.  The  axial  strain  e  way  be  expressed 
in  terms  of  the  strains  in  ( 14)- ( 16)  using  the  second-order  tensor 
transformation  rule, 

2  2 

£x  =  cos  6  ej  +  sin  e  -  sin6  cose  £g  (18) 

Substitution  of  (14)-(17)  into  (18)  yields 

S  ’  cx  "  0x/oo  (19) 

where  e®  is  the  elastic  axial  strain,  and 

4  2  2  !- 

h  i  (sin  e  +  a^g  sin  8  cos  e)z  (20) 

Observe  also  from  (13)  and  (17)  that 

°o  =  hox  <21> 

We  can  obtain  the  function  h(e)  used  by  Sun  and  Chen  by  multiplying  (20)  by 
/3/2.  Equation  (19)  is  the  same  as  derived  by  them  from  a  plasticity 
model  for  loading  behavior;  this  model  will  be  discussed  in  Section  3. 

Experimental  information  on  c  -  ox  behavior  for  two  fiber  angles  e  may  be 
used  with  (19)  to  evaluate  x  a^g  and  the  function  e  =  e  (o  )  . 
(Alternatively,  one  may  use  data  from  several  fiber  angles  to  determT ne°the 
a66  which  minimizes  the  data  spread  in  the  e  (o  )  plot.)  Results  from 
tests  at  other  fiber  angles  then  serve  to  check  (ly).  A  simple  power  law 

co  =  A  °on  <22> 

where  A  and  n  are  positive  constants,  was  reported  by  Sun  and  Chen  to  fit 
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the  data  out  to  specimen  failure  (e  -  1%);  for  a  boron/aluminum  composite 
n  =  5.8  and  agg  =  4,  whereas  for  gr$)hite/epoxy  n  =  3.7  and  agg  =  2.5.  (The 
constant  agg  used  by  Sun  and  Chen  is  one-half  of  the  agg  used  here.) 
Values  of  agg  =  4  and  n  =  2.4  have  been  obtained  recently  Dy  Mignery  and 
Schapery  (1988)  from  studies  of  unidi rect ional  and  angle-ply  laminates  of 
the  same  rubber-toughened  graphite/epoxy  material  used  to  develop  the 
curves  in  Fig.  1.  Although  the  latter  exponent  (n  =  2.4)  is  smaller  than 
that  reported  by  Sun  and  Chen  for  an  untoughened  unidirectional 
graphite/epoxy  material  (n  =  3.7),  the  angle-ply  stress-strain  curves 
(Mignery  and  Schapery,  1988)  exhibit  a  larger  degree  of  nonlinearity 
because  the  total  axial  strain  range  is  approximately  5x,  as  compared  to  1% 
in  the  former  study. 


In  the  much  earlier  work  of  Lou  and  Schapery  (1971)  it  was  found  that  the 
parameter  o  in  (13)  accounted  for  the  effect  of  stress  state  on  the 
functions  us°eo  to  characterize  nonlinear  viscoelastic  behavior  of  a 
glass/epoxy  composite.  The  motivation  for  the  use  of  this  parameter  came 
in  part  from  the  observation  that  the  octahedral  shear  stress  t  .  can 
normally  be  used  to  correlate  multiaxial  yielding  of  plastics  (just  as  for 
metals).  As  a  simplification,  the  matrix  was  viewed  as  a  uniformly 
stressed  layer  of  material  sandwiched  between  layers  of  rigid  fiber 
material;  i.e.,  the  lines  in  Fig.  2  at  the  angle  0  were  imagined  to  define 
layers  rather  than  fibers.  Using  the  principal  material  axes,  Fig.  2,  this 
shear  stress  is 


‘oct 


(°i“ 


>)2+ 


-~33)‘ 


[a-. 


6(°42+ 


+  0 


-  2)}H 


(23) 


where  the  in  this  equation  are  the  stresses  in  a  matrix  layer. 

For  a  matrix  in  plane  stress  o?  and  o ,  are  the  same  as  the  stresses 
o 2  and  or  acting  on  a  composite  consisting  of  parallel  layers  of  matrix 
ano  reinforcement  material.  A  factor  v  was  also  introduced,  as  defined 
by  the  relationship  o,  =  v  .  For  a  linear  elastic,  isotropic  matrix 
v  is  the  Poisson's  ratio,  and  for  an  incompressible  elastic  or  rigid- 
plastic  matrix  vg  =  0.5.  Use  of  these  idealizations  in  (23)  yields 

Toct  =  (2/3c)iS  (o22+  Co62)H  (24a) 

where 

c  =  3/(1  -  ve  +  ve2)  (24b) 


As  reported  by  Lou  and  Schapery  (1971)  a  finite  element  analysis  of  a 
linear  elastic  composite  with  a  square  array  of  fibers  was  made  to  predict 
the  average  octahedral  shear  stress  in  the  matrix.  Apart  from  a  numerical 
factor,  (24)  was  found  to  be  a  fairly  good  approximation  to  this  average. 
Considering  c  to  be  the  arbitrary  constant  agg,  it  is  seen  that  (24a)  and 
(13)  are  equivalent  parameters  for  characterizing  nonlinear  behavior.  It 
is  also  of  interest  to  find  from  (24b)  that  c  =  4  when  m  =  h  and  c  =  3.88 
when  v  =  0.35;  the  former  value  is  the  same  as  found  Experimental ly  by 
Sun  ande  Chen  (1987)  for  the  boron/aluminum  composite  and  by  Mignery  and 
Schapery  (1988)  for  the  rubber-toughened  graphite/epoxy  composite;  the 
latter  value  of  c  was  reported  by  Lou  and  Schapery  (1971)  for  glass/epoxy 
material . 


Most  of  the  experimental  work  reported  above  is  for  proportional  loading, 
(17).  However,  that  of  Mignery  and  Schapery  (1988)  involves 
nonproportional  loading  of  the  plies  in  an  angle-ply  layup.  Ti.^e  studi>_o 
provide  limited  experimental  support  for  (11).  We  are  currently  making 


8 


additional  studies  of  angle-ply  and  unidirectional  laminates  under  loading, 
unloading  and  reloading  to  address  the  applicability  of  (10)  and  (42)  for 
toughened  and  untoughened  graphite/epoxy  composites. 

It  should  be  observed  that  the  difference  between  loading  and  unloading 
curves  in  the  model  (6)  is  characterized  by  one  scalar  factor 
dP/do  ,  where  P  =  P(o  ,  S).  The  loading  curves,  do  /dt>0  .are  predicted  by 
using  o  =  o  .  For  unloading,  do  /dt<0  ,  the  thermodynamic  requirement  of 
positive  entropy  production  and  one  path-independence  of  the  unloading  work 
are  violated  unless  S  is  constant  Schapery  (1988).  Consequently,  fur 
arbitrary  stress  histories,  S  is  always  the  largest  value  of  a  up  to  the 
current  time. 

This  representation  does  not  account  for  the  difference  between  unloading 
and  reloading  curves.  Tonda  and  Schapery  (1987)  were  able  to  account  for 
this  difference  for  an  untoughened  graphite/epoxy  composites  using  linear 
viscoelasticity  theory;  the  approach  to  combining  the  effects  of 
viscoelasticity  and  structure  changes  was  developed  earlier  (Schapery, 
1981).  Whether  or  not  this  approach  is  able  to  account  for  all  of  the 
hysteresis  is  not  presently  known.  It  may  be  necessary  to  introduce 
another  S-parameter  which  is  activated  at  the  start  of  reioaoing. 

3.  THE  NORMALITY  RULE  FOR  INELASTIC  STRAINS 


Let  us  now  compare  the  normality  rule  employed  in  plasticity  theory  to 
predict  plastic  strain  increments  with  the  type  of  normality  contained  in 
(4).  Following  Sun  and  Chen  (1987),  we  take  =  k  as  the  yield  condition, 
k  ic  a  scalar  that  varies  with  the  mount  of  plastic  straining. 
Plastic  strains  are  introduced  in  the  same  way  as  is  commonly  done  for 
metals, 

5  d=i  - <«> 

where  deP,  de.,  and  de*:  are  infinitesma!  changes  in  plastic,  total,  and 
elastic  strains,  respectively.  The  elastic  strains  are  assumed  to  be 
linear  in  the  stresses. 


c>  =  Vj 


(26) 


where  are  the  constant  compliances.  The  associated  flow  rule  for 

plastic  strain  increments  is 


r,  3°o 

de?  =  dx 

1  3o.j 


(27) 


where  dx  is  a  scalar.  This  equation  shows  that  de?  is  a  vector  which  is 
normal  to  the  surface  aQ  =  constant.  From  (8)  and  (27), 


d£P  -  2ai jdjdx 


(28) 


F  o 
be 


r  proportional  stressing  o.  =  k^o  (where  the  k^  are  constants)  (28)  may 
integrated  to  obtain?the  total  prastic  strains, 


n  3°n 

e?  =  7 - ( f o  d x ) /  o 

1  9o.j  v  0  '  0 


(29) 


which  is  also  a  vector  normal  to  the  surface  o  =  constant.  The  total 

strain  is 


e 


i 


P 

i 


(30) 
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which  may  be  compared  to  the  strain  (6)  derived  from  a  dual  strain  energy 
density.  The  "inelastic  strain"  vector  in  (6), 


I 


3P_ 

3  o , 


3o, 


(31) 


is  normal  to 
contrast  to 
proportional 


1  '  3oo  3oi 

surface  o 

t.,  the  normality  of  el  exists  for  proportional  and  non¬ 
stressing.  Observe  also  that  this  normality  is  preserved 
during  unloading  and  reloading;  recall  that  the  coefficient  aP/3o  depends 
on  both  o  and  S,  and  that  S  =  o  only  when  o  is  equal  to  its  largest 
value  (considering  all  values  up  to  The  current  time). 


3  =  constan^,  just  as  in  (29).  However,  in 
iTity  of  e )  exists  for  proportional 
also  that  this 


Consider  next  for  further  comparison  a  type  of  normality  discussed  by  Rice 
(1971)  for  incremental  inelastic  strains.  He  developed  (4)  from 
thermodynamics  with  internal  variables  and  used  it  in  a  study  of  inelastic 
behavior  of  metals;  S  is  one  of  possibly  many  internal  variables.  A  change 
in  strain  due  to  infinitesimal  changes  in  both  a.  and  S  is,  from  (4), 


where 


3  0 

3o^30j^°j 


sG 


3  0  . 


-dS 


(32) 


G  =  - 3w  /aS  (33) 

o 

Rice  observed  that  when  elastic  and  inelastic  strains  are  defined  through 
increments,  as  expressed  by  the  first  and  second  terms  in  (32), 
respectively,  the  incremental  inelastic  strain, 

cki  s  -ff^S  (34) 

is  normal  to  the"yield"  surface  G  =  constant.  In  fracture  mechanics  G  (33) 
is  called  the  "energy  release  rate".  When  there  are  two  or  more  structure 


parameters 

S  (m  = 
m 

-{■ 

1,  2,  . . . ) , 

3G 

where 

—nLdSm 

3o^  m 

(35) 

-3W  / 3 S_ 
a  m 

(36) 

Thus,  the  mth  component  of  dej  is  normal  to  the  respective  surface,  Sm  = 
constant,  as  noted  by  Rice. 


When  we  use  the  special  form  for 
inelastic  strains  become  ~ 

3  w 


w  in  (5), 

o 


Rice's  incremental  elastic  and 


32w 


aO 


3o  .  3a  • 

*  J 


sG 


3  o 


dS  =  - 


J 

32P 


3  o  .  3  o  • 
'  J 


-da . 


32P 


the 


3  S  3  c 


-dS  =  - 


32P 


3  o  ^  3  a  j 


-d  a  . 


(37) 


3  a. 


■ I 


3  S3  c 


dS  jf  ^ 


Notice  that 
increment  in 
term  in  (37). 
matrix  -3w  /da. do 


is  normal  to  the  surface  o  =  constant  and 
elastic  strain  defined  in  (7)  ts  equal  to  only 
Observe  also  that  the  tangent  elastic 
used  in  defining  the  incremental  elastic 


that  an 
the  first 
compl i ance 

...  , _  _  ...  __  v  ....  _ _  _ _  strains  in 

(32)  is  a^unc’tic^i  of  the  structure  parameter  S  as  well  as  stresses,  while 

that  based  on  the  elastic  strain  in  (7),  -a  -  J - “-1-  --  iU~ 

stresses. 


w  /3o.3o.,  depends  only  on  the 
oo  i  J 
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4.  A  CONSTITUTIVE  EQUATION  WITH  LINEAR  UNLOADING  BEHAVIOR 

In  characterizing  the  effect  of  damage  on  composite  material  behavior,  it 
is  commonly  assumed  that  the  material  is  linearly  elastic  when  damage  is 
constant.  This  linearity  assumption  is  equivalent  to  using  a  dual  strain 
energy  density  in  which  stress  dependence  is  limited  to  first  and  second 
order  terms. 


w 

0 


-b. 


Vi 


bijVj 


(39) 


where  bQ,  b^  and  bjj  may  be  functions  of  one  or  more  structure  parameters 
$m.  In  this  case  t*4  strains  (4)  are 


=  bi 


b  .  .0  . 
ij  J 


(40) 


The  residual  strains  and  compliances  bi  j  may  vary  with  stress  history 
through  changes  in  Sm;  only  one  S  will  beJused  here.  The  strain  energy 
density  is  related  to  w^  through  (3), and  may  be  written  as 


w 


co  +  ci 


(41) 


which  provides  the  stresses 


°i 


=  c.  +  c. ,  c,  (42) 

be  obtained  by 
be  sufficiently 

stress-strain  behavior  if  the  unloading  and 
approximated  by  the  same  straight  line  whose 
S  (as  shown  in  Fig .  3) . 


iJ  J 

The  relationship  between  the  b's  and  c's  may  of  course 
comparing  (40)  and  (42).  These  second-order  energies  may 
general  to  predict  ply 
reloading  curves  can  be 
position  (c.j)  and  slope  (c^)  vary  with 


The  work  (Jo.de.)  and  dual 
n^rocrgs  a  re  1  i  nd  ep  ende  nt  of 


proces 

1988), 


work 
path  or 


(-  /e-do,)  during 
history  'if  and 


3w_ 

IT 


or 


3w 


sS 


=  9 


structure-change 
only  if  (Schapery, 


(43) 


Fig.  3.  Stress-strain  behavior  according  to  (42),  showing 
loading,  unloading  and  reloading. 


where  g  is  at  most  a  function  of  S;  the  quantity  g  is  the  specific  fracture 
energy  if  S  is  the  fracture  surface  area  of  a  crack.  (Equation  (43)  is  not 
limited  to  the  second-order  energies  (39)  and  (41).)  As  shown  by  Schapery 
(1988)  S  can  always  be  chosen  so  that  (43)  reduces  to 


Observe  that  the  term 
g  in  (43). 


(=  -  b  )  can  be  omitted  as  it  can  be  absorbed  in 


It  should  be  added  that  the  derivatives  aw/aS  and  aw/aS  are  always  equal 
which  may  be  easily  shown  by  taking  the  differential  of  (3).  Equation  (44) 
provides  the  relationship  for  predicting  S  as  a  function  of  stress  or 
strain.  Thermodynamic  theory  requires  dS/dt  >  0  (Schapery,  1988);  thus,  if 
(44)  predicts  dS/dt  <  0,  S  is  actually  constant  and  (44)  is  to  be 
di  sregarded. 

For  the  second  order  energy  (39)  with  bQ  =  0,  the  equation  for  S  is 


-vr  1 


Although  (39)  is  only  of  second  order  in  the  stresses,  it  is  still 
sufficiently  general  to  mathematically  represent  Sun  anj  Chen's  data 
discussed  in  Section  2.  Indeed,  this  may  be  done  by  assuming  the  b^  are 
constants  and  then  using 


BSaij 


where 


n  -  1 
n  +  1 


B  «  KUr(2,r)r 


Also,  S^  s  are  the  c 
constants  appearing 
(45)  and  (46)  yield 


are  the  constant  elastic  compliances. 


(22);  observe  that 


(46a) 


(46b) 

n  are  the 
Equat ions 


=  (Br/2 ) 


(n+1) /2  (n+1 ) 


During  loading,  do  /dt  >  0,  (47)  is  used  in  (46a)  to  predict  instantaneous 
values  of  b^.  Tor  unloading,  do  /dt  <  0,  the  coefficients  b^  are 

constant  because  S  has  a  constant  va?ue  equal  to  that  at  the  start  of 
unloading.  Upon  reloading,  S  again  changes  in  accordance  with  (47)  when 
a  reaches  its  largest  past  value.  Unloading  and  reloading  data  are  not 
reported  by  Sun  and  Chen  (1987),  and  thus  the  range  of  applicability  of 
this  particular  model  cannot  be  assessed  at  this  time.  It  is  important  to 
notice  that  this  phenomenological  characterization  is  not  necessarily 
limited  to  brittle  or  to  ductile  composites,  as  Sun  and  Chen's  results  are 
for  both  types. 

Finally,  we  should  mention  that  the  theory  based  on  path-independence  of 
work  has  been  successfully  employed  in  limited  studies  of  particle- 
reinforced  rubber  (Schapery,  1987b),  and  a  thermoplastic  composite  (Dan 
Jumbo  et  al . ,  1987).  In  the  former  case  nearly  all  nonlinear  behavior  was 
expressed  in  terms  of  S-dependence  of  b ^  ;  in  the  latter  case  the  residual 
strains  b^,  instead  of  b^j,  were  used  to  account  for  most  nonl ’nearities. 
The  small  amount  of  nonlinearity  that  was  not  adequately  represented  by  the 
second-order  energy  functions  was  apparently  due  to  the  large  strains 
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(  =  60^)  in  the  filled  rubber  specimens  and  fiber  or  microfibril  alignment 
(causing  an  increase  in  modulus  for  loading  in  the  fiber  direction)  in  the 
thermoplastic  composite. 

5.  CONCLUDING  REMARKS 

A  possible  approach  to  predicting  multidirectional-fiber  laminate  behavior 
would  consist  of  using  a  unidirectional  ply  energy  density,  such  as  given 
by  (5)  or  (41),  with  the  usual  displacement  assumptions  of  lamination 
theory  (Christensen,  1979).  Delaminations  and  their  growth  could  be 
accounted  for  essentially  in  the  same  way  as  done  for  linear  and  nonlinear 
elastic  laminates,  but  with  additional  bookkeeping  when  there  is  any 
appreciable  difference  between  loading  and  unloading  stress-strain 
behavior.  The  work  of  deformation  (which  is  equal  to  w.  ;  w  +  S  if  the 
second  equation  in  (44)  is  used  to  predict  S)  is  treated  just  like  strain 
energy  in  nonlinear  elastic  fracture  mechanics  (Schapery , 1987a) ;  in 
particular,  w-.  is  used  in  strain  energy  release  rate  and  J  integral 
calculations. 

With  brittle-matrix  composites,  a  significant  number  of  transverse  ply- 
level  cracks  may  develop  prior  to  structural  failure  (Johnston,  1987). 
These  cracks  are  somewhat  planar  with  the  plane  parallel  to  the  fibers  and 
perpendicular  to  the  lamination  plane.  Typically,  after  rapid  growth,  they 
are  arrested  at  the  ply  boundaries.  If  more  than  one  fiber  orientation  is 
used,  a  laminate  usually  is  capable  of  supporting  loads  well  above  that  at 
crack  initiation.  Whether  or  not  one  S-parameter  is  sufficient  to  account 
for  a  general  type  of  inelasticity  which  includes  transverse  cracks 
requires  further  study.  It  should  be  observed  that  even  with  only  one 
parameter,  an  appreciable  effect  of  these  cracks  on  the  laminate  behavior 
may  be  taken  into  account  through  the  way  w  or  w^  depends  on  S;  for 
example,  by  may  have  the  form  in  (46a)  at  small  S,  and1  then  a  considerably 
rliffprent  Hrm  at  large  S  when  transverse  cracks  develop.  Physically,  S 
may  reflect  micro-damage  (e.g.  rubber  particle  cavitation)  and  plastic 
deformation  until  transverse  cracks  develop,  and  then  at  larger  S-values 
account  for  these  mechanisms  as  well  as  transverse  crack  density.  If  the 
effect^  of  crack  density  and  its  growth  are  not  sensitive  to  properties  of 
adjacent  plies  with  different  fiber  angles,  an  experimental  program  could 
use  the  simple  angle-ply  layup.  Similar  observations  can  be  made  for 
distributed  interior  delaminations  (Harris  et  al . ,  1987);  however,  at  least 
two  plies  would  comprise  the  basic  element  of  a  laminate. 

We  are  presently  using  these  ideas  to  characterize  and  predict  the 
mechanical  response  of  untoughened  and  toughened  graphite/epoxy  laminates, 
recognizing  that  the  proposed  method  has  to  be  considered  as  tentative 
until  a  significant  amount  of  additional  experimental  and  analytical 
studies  are  made.  Such  studies  should  help  to  establish  the  range  of 
validity  of  the  work-potential  method  as  well  as  define  the  experimental 
program  needed  for  a  complete  characterization.  Micromechanical  models  of 
damage  in  linear  elastic  composites  (Wang  and  Haritos,  1987)  should  be 
helpful  in  analytically  modeling  the  effect  of  distributions  of  cracks  on 
moduli  or  compliances,  and  thus  reduce  the  experimental  effort.  Schapery 
(1987b)  used  this  approach  in  an  elementary  model  to  relate  the  orthotropic 
elastic  properties  of  a  particulate  composite  to  a  statistical  distribution 
function  which  characterized  the  damage,  and  employed  an  evolution  equation 
like  (44)  to  predict  the  change  in  properties  through  an  S-parameter  which 
is  an  overall  measure  of  the  damage.  A  similar  procedure  should  be 
applicable  to  laminates. 
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ABSTRACT 

A  method  of  laminate  characterization  and  analysis 
is  descriDed  in  which  growing  damage  and  other  inelastic 
phenomena  are  treated  using  the  same  mathematical 
formalism,  thus  simplifying  the  description  of 
mechanical  response.  It  is  based  on  the  observation 
that  the  applied  work  is  not  sensitive  to  many  details 
of  the  deformation  history.  Following  a  brief 
discussion  of  the  thermodynamically-based  theory,  we  use 
a  special  version  along  with  experimental  data  on  a 
graphite/epoxy  composite  to  obtain  an  explicit 
mathematical  characterization  of  a  unidirectional  ply. 
Predictions  of  mechanical  response  are  then  compared  to 
experimental  results  for  a  variety  of  layups,  one  of 
which  delaminates  from  the  edges. 
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Delamination  depth 
Coefficients  in  strain  energy 
Young's  Modulus  GPa  ( tb/in^xlO^) 
Thermodynamic  force 
Strain-dependent  coefficients 
Shear  modulus  GPa  (lb/1n^xl0®) 

Reduced  modulus 

Structural  parameter 

Work  needed  for  structural  changes/vol 

Cube  root  of  S 

Strain  energy/vol 

Total  work  input/vol 

Total  work  input 

Cartesian  coordinate 

St ra  in 

Poisson' s  ratio 
Stress  (lb/1n?xl03) 


1.  INTRODUCTION 

The  nonlinear  behavior  of  unidirectional  and 
multidirectional  fiber  composites  traditionally  has  been 
modeled  using  elasticity  theory  (1)  or  plasticity  theory 
in  which  the  unloading  moduli  are  constant  (2,3). 
iiowever,  resin  matrix  composites  often  exhibit  nonlinear 
behavior  which  is  due  at  least  in-part  to  inelastic 
mechanisms  that  alter  the  unloading  moduli.  As 
1 1  lustrated  in  F ig .  1  for 


Fig.  1  Stress-strain  behavior  of  Hercules  AS4/3502 
angle-ply  laminate  |t30|j^  showing  inelastic 
behavior. 


the  brittle-resin  composite  Studied  later  in  this  paper, 
the  inelasticity  (i.e.  strain  history-dependence)  is 
neither  small  enough  to  neglect  nor  large  enough  to  use 
classical  plasticity  theories  in  which  unloading  follows 
the  initial  modulus  EQ.  For  several  unidirect  iona  1  and 
angle-ply  layups  of  this  graphite/epoxy  material  we  have 
found  that  the  residual  strain  upon  load  removal  is 
typically  20-40%  of  that  for  unloading  along  the  E0- 
line.  There  is  also  a  small  amount  of  rate  or  time- 
dependence  even  in  the  room  environment. 

In  this  paper  we  discuss  a  way  of  characterizing 
nonlinear  inelastic  behavior  that  may  arise  from  a 
variety  of  mechanisms,  including  microcracking, 
delamination,  void  growth,  shear  yielding,  and 
crystalline  slip,  and  therefore  is  not  limited  to  resin- 
matrix  or  metal-matrix  composites.  It  is  based  on  the 
theory  in  (4),  which  uses  Rice's  (5)  thermodynamic 
description  of  inelastic  behavior.  Some  motivation  for 
the  theory  has  come  from  our  experimental  studies  of 
rubber-toughened  and  untoughened  graphite/epoxy  (6,7) 
under  axial -torsional  loading,  which  show  that  the 
stresses  and  mechanical  work  are  practically  independent 
of  deformation  history  for  suitably  limited  paths.  This 
limited  path- independence  leads  to  a  mathematical 
description  of  mechanical  behavior  which  is  analogous  to 
that  used  for  predicting  stable  crack  growth  and  its 
effect  on  global  structural  response. 

In  Section  2  the  theory  for  a  unidirectional  ply  or 
laminate  is  outlined,  and  the  plane-stress  case  used  ’n 
later  sections  is  described;  the  theory  allows  for 
elastic  nonlinearity  (such  as  that  due  to  fiber- 
straightening)  and  inelastic  nonlinearity  during 
loading.  Sections  3-5  discuss  the  experimental  program 
and  both  linear  and  nonlinear  behavior;  after  the  basic 
ply  characterization  is  accomplished,  response  of 
several  laminates  is  predicted  and  shown  to  be  in  good 
agreement  with  the  experimental  findings. 
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2.  CONSTITUTIVE  EQUATIONS  FOR  A  UNIDIRECTIONAL  PLY 
Basic  theory 

The  coordinate  notation  for  a  unidirectional 
laminate  or  ply  is  shown  in  Fig.  2;  the  x.  axis  is 
parallel  to  the  fibers  while  the  *3  axis  is  normal  to 
the  ply  plane.  The  stresses  Oj  and  strains 


Y 


t,  (1*1,2, ...6)  are  referred  to  the  principal  material 
coordinates  x^(k  =  l  ,2,3).  We  also  employ  the  standard 
notation  in  which  i=4,5,6  refer  to  the  shearing  stresses 
and  strains,  where  these  shearing  strains  are  twice  the 
tensor  components. 

The  stresses  are  related  to  strains  through  the 
strain  energy  density  (per  unit  Initial  volume)  w  = 
w(cf.  SJ. 


°i  =  3w^ci  (1) 

where  the  Sm  are  so-called  "structural  parameters". 
Temperature  and  other  parameters  (such  as  moisture)  may 
enter;  but  for  simplicity  we  assume  they  are  constant, 
and  thus  do  not  explicitly  show  them  as  arguments.  As 
n'?ny  sm  4S  needed  are  used  to  account  for  the 
mi crostructura  1  changes  which  produce  the  Inelastic 
behavior. 

It  will  be  helpful  to  also  use  the  dual  strain 
energy  density  w  *  w  (o.,  S  ) , 
a  oi  m' 9 


w0  =  w  -  °i  £i  (2) 

which  Is  seen  to  be  the  negative  of  the  complementary 
strain  energy  density.  The  summation  convention  is 
employed  wherein  a  repeated  Index  implies  summation  over 
its  range.  By  using  Eq.  (1)  with  the  differential  of 
Eq.  (2)  for  constant  Sm,  we  find 


C1  1  -3V3o1  (3) 

Changes  in  S_  may  be  related  to  stresses  or  strains 
by  using  a  constitutive  equation  for  the  thermodynamic 
forces  fm;  these  forces  are,  by  definition, 

fm  5  -3w/3Sm  <«) 

The  differential  of  Eq.  (2),  for  independent  changes 
dSm,  yields 
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The  second  law  of  thermodynamics  allows  only  those 
changes  for  which 
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where  the  overdot  denotes  a  derivative  with  respect  to 
time.  As  the  constitutive  equation  for  f  we  specify 

that  for  each  active  parameter,  i.e.  5  »  0? 
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3ws/35m 


(7) 


where  ws  =  ws(Sm)  is  a  constitutive  function  of  S. ; 

also,  for  those  m  in  which  f  «  aw  /sS  then  5  =  0. 

Each  fm  is  viewed  as  the  fofte  avli  la'Ele  to  Produce 

changes  in  the  associated  Sm,  while  aw  / aS  is  the  force 
required  for  these  changes.  The  sofijtiPli  of  Eq.  (?) 

yields  S  =  S  or  S  «  S  (a,)  for  each  active 

parameter.  m  m  i 

The  total  work  input  per  unit  initial  volume  during 
actual  elastic  or  inelastic  processes  is 


It  is  always  possible  to  select  c.  and  Sm  such  that  they 
vanish  in  the  initial  or  reference  state,  and  therefore 
we  use  such  a  choice  throughout  this  paper.  Given  Eq. 
(7),  one  may  easily  show  that  wT  is  a  potential  function 
of  the  state  (tj(  S^)  and,  in  particular,  that 


WT  *  w  ♦  -s  .  V  3wT(cj'Sm(tj^/3cl  (9> 

(Without  loss  in  generality  we  specify  that  w  =  w  =  0 
in  the  initia1  state.)  According  to  Eq.  (9),  the  total 
work  consists  of  the  work  of  straining  w  plus  the  work 
of  structural  change  ws .  Moreover,  from  Eq s . ( 6 )  and 


i  0  (10) 

Quadratic  strain  energy  functions 

Consider  now  as  a  candidate  dual  strain  energy  the 
quadratic  form 

wo  *  -v  Vi  -!bijVj  +  ?v(V  V(sr  “jH11) 

The  last  term  allows  for  up  to  six  independent 
structural  parameters,  while  all  coefficients  may  depend 
on  these  as  well  as  other  parameters;  unless  this 
dependency  is  restricted  in  some  way  the  last  term  is 
redundant.  The  strain  Eq.  (3)  becomes 
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(12) 


Both  b^  and  r..  contribute  to  the  residual  strains  (i.e. 
the  strains  when  o.  =  0).  As  a  further  specialization 
of  the  theory,  we  assume  that  the  bj  are  independent  of 
Sm  and  and  that  the  remaining  coefficients  b.,.  and 
r^j  depend  on  only  one  parameter,  S,  which  we  take  io  be 
w  .  The  term  bQ  affects  only  Eq.  (7),  and  can  be 
absorbed  Into  w$,  although  the  inequality  (10)  would 
still  apply  only  to  w$ .  When  the  Sm  define  only  the 
state  of  damage  in  the  form  of  cracks  and  voids,  b  is 
the  surface  free  energy,  and  thus  may  be  viewed  as  the 
energy  available  for  healing;  here  we  assume  b  =  0  so 
that  ‘  >  0.  0 

With  these  restrictions  on  the  coefficients,  Eqs. 
(5)  and  (7)  yield 

S  *  0  (13) 


and 
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rij(Sj-  Oj)  =0  if  s,  -  0  (14) 

If  all  *  0  and  if  r^  is  positive  definite,  then  Eq. 
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(14)  implies 


SJ  " 

and  Eq.  (12)  reduces  to 

‘i  *  bi  +  b< 
while  Eq.  (13)  becomes 
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We  may  now  give  a  physical  interpretation  for  the 
three  sets  of  coefficients  in  Eqs.(ll)  and  (12).  The 
b;  ;  are  compliances  that  in  general  vary  with  the  work 
done  in  changing  the  microstructure,  S.  Since  S=  ws, 
Eq.  (10)  implies  5  >  0.  If  we  know  the  functions 
b,  i( S) .  then  Eq.  (17)  may  be  used  to  predict  S  in  terms 
ofJstress  history;  if  the  equation  predicts  S  <  0,  then 
we  set  S  «  S  ,  where  S  is  the  largest  value  of  S  up  to 
the  current 8  time.  1  If  Eq.  (17)  is  subsequently 
satisfied,  then  again  it  is  used  to  predict  S. 

Whenever  S.  •  a.,,  then  instead  of  Eq.  (16)  we  have 
from  Eq.  (12),  J  J 


■i  =  v  rijsj  *  (V  rij)3j 


(18) 


It  is  seen  that  by  means  of  r^,  we  may  simulate  the 
effect  of  internal  surface  roughness  and  other 
irregularities  In  resisting  crack  closing  and  sliding 
and  void  collapse  during  unloading;  whenever  all 
S.  =  o  .  these  coefficients  have  no  effect  on  the  strain, 
Eq.  (lt>).  Had  we  allowed  for  dependence  of  on  Sm, 
then  these  coefficients  would  produce  the  type  of 
residual  strains  usually  associated  with  plastic 

deformation  mechanisms,  such  as  crystalline  slip. 

The  strain  energy  density  w  corresponding  to  Eq. 
(11)  may  be  obtained  from  Eq.  (2)  after  stresses  are 
expressed  in  terms  of  strains.  Whether  or  not  the  r^; 
enter,  it  is  of  the  form 
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From  Eq.  ( 1 )  , 
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By  comparing  Eq.  (20)  to  (16)  or  (18),  we  may  obtain  the 
coefficients  in  Eq.  (20)  in  terms  of  those  In  the 
former  equations.  Figure  3  illustrates  the  loading 


and  unloading  behavior  in  terms  of  the  coefficients  in 
Eq.  (20),  where  S  represents  all  structural  parameters. 

In  all  of  the  subsequent  work  we  shall  be  concerned 
with  only  loading  processes,  defined  here  as  those  for 


which  $  >  0  as  well  as  S 
Moreover,  we  shall  assume  b 
discussion  on  thermal 
simplifications 
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or  else 


.  =  o,  or  eiae  r^j 

■  —  c i  —  0 ,  except  inJa  brief 
residual  strains.  With  these 
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Also  noting  that  Eqs.  (4) 
3w/aS  =  aw  /aS,  then  from  Eq.  (13), 
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From  Eq.  (9) 


=  -1 


(21) 

(22) 

(5)  imply 

(23) 


and 


S  *  Wj  -  w  (24) 

which  shows  that  the  Instantaneous  value  of  S  is  the 
shaded  area  in  Fig.  3. 


Constitutive  equations  for  plane  stress 

Constitutive  Eq.  (22)  has"  the  same  form  as  that  for 
a  linear  elastic  material.  If  we  assume  the  material  is 
orthotropic,  with  principal  material  directions  xk,  then 
we  may  use  standard  linear  elastic  equations  to 


characterize  ply  behavior 
moduli.  From  (8), 

in  terms  of 

principal 
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(26a) 

°12  =  V1 2Q22 

• 

D  =  1  -  v12v?1 

(26b) 

V21  =  V1 2E2/E  1  ’ 

Q66  *  G12 

(26c) 

Only  the  Young's  modulus  for  loading  normal  to  the 
fibers  Ep  and  the  shear  modulus  Gjp  are  assumed  to  vary 
as  a  result  of  inelastic  mechanisms.  The  remaining  two 
principal  properties,  i.e.  Young's  modulus  in  the  fiber 
direction  E,  and  Poisson’s  ratio  v.,,  are  assumed  to  be 
independent  of  S,  which  is  reasonable  if  the  fibers  do 
not  break.  Indeed,  our  experimental  results  from 
loading  and  unloading  tests  parallel  to  the  fibers  sho~ 
negligible  hysteresis  in  axial  and  transverse  strains, 
even  when  the  maximum  strain  is  close  to  the  ultimate 
value.  However,  there  is  some  elastic  nonlinearity 
which  is  no_t  negligible  for  our  purposes;  Ej  increases 
and  Eo  and  v,?  decrease  with  increasing  strain. 

In  order  to  account  simultaneously  for  inelastic 
effects  in  Ep  and  Gj^  and  nonlinear  elastic  effects  in 
Ej,  Ep  and  v12  we  propose  the  following  modified  form  of 
Eq.  (25), 
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Q 1  j  E  j  c  j  *  v<3?2f  12  ^ 2 
vQ22f 21  1  *  Q22f2c2 


(27a) 

(27b) 
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where 

1 12 

=  G12'12 

(27c) 

V  = 

and 

v12(°). 

*  1  =  /  E12£*i:1  *  *2  ;  ^  E2dc2 

0  0 

(28) 

fi  * 

Eo^l> 

f  1  dt2  r  E90(c2) 

(29) 

eTToT 

’  12  -  v  dcj  ’  2  -  E2(0) 

The  strain-dependent  functions  fj  and  f^  are  the  ratios 
of  the  secant  modulus  E  -  o./t.  and  tangent  Poisson's 
ratio  at  a  given  axial  strain  to  those  at 
S-Cj=0(for  o,  *  0,  Op  «  T..  =  0).  Allowance  has 
been  made  for  elastic  rfbnl  ineari  ty  under  transverse 
loading  (o,  <  0,  o.  «  t._  =  0)  through  f 2;  the  quantity 
Eqg  in  Eq.  (29)  is  the  Secant  modulus  for  S  =  0  from  the 
90  specimen.  Observe  that  Ij  =  -  c-(t.)/v,  where  this 
<-2  is  the  transverse  strain  in  a  uniaxial  test  with 
o.  <  0,  o,  -  i 12  =  0- 

The  r actor  D  in  Eq.  (26b)  is  practically  unity,  and 
although  it  is  not  constant  we  may  use  its  initial 
linear  elastic  value  in  all  work  without  any  significant 
error.  All  strain  dependence  of  the  moduli  Is  now 
accounted  for  through  fj,  f 2,  and  fio,  and  therefore  Qii 
and  Q22  in  Eq.  (27)  are  those  in  Eq.  (26a);  vi2.,  £1  and 
Qjl  are  constants,  while  Q 22  varies  with  S  through  E2; 
E 2  and  G 1 o  n0w  vary  only  with  S. 

The  form  of  Eq.'  (27a)  and  Eq.  (27b)  is  such  that 


which  is  a  necessary  and  sufficient  condition  for 
construction  of  the  strain  energy  function  in  Eq.  (1). 
We  find  that 

w  *  QuIn  +  $22l22 

+  vQ22!2ii  +  G12> 12/2  <31) 

where 

*11  "  ^  clEldtl  ’  *22  5  ^  c 2 f 2c*t 2  ^33^ 


supplied  In  unidirectional  prepreg  form  and  cured  in  our 
air-cavity  press  using  the  supplier's  specification  for 
the  cure  cycle.  The  specimens  were  cut  from  12"xl2" 
plates  and  stored  in  a  dessicant  until  they  were  strain- 
gaged  and  tested  under  uniaxial  loading  in  an  amDient 
environment  of  about  7 5 * F  and  at  a  tensile  strain  rate 
of  0.005/min. 

All  unidirectional  and  angle-ply  layups  consisted 
of  12  plies.  Most  specimens  were  0.5  inch  wide,  with  a 
length  of  7.5  Inches  between  the  glass-epoxy  end-tabs. 
Some  1x10  inch  specimens  were  also  tested  for  comparison 
purposes;  no  significant  differences  between  the  0.5  and 
1.0  inch  wide  specimens  were  found.  By  screening  out 
all  samples  whose  thicknesses  differed  by  more  than  5 /., 
very  little  specimen-to-specimen  differences  in 
mechanical  behavior  were  observed.  In  most  cases  the 
results  presented  are  for  an  average  of  two  specimens. 
The  off-axis,  unidirectional  specimens  exhibited  the 
greatest  variability,  but  we  have  not  attempted  to 
quantify  the  scatter.  Two  pairs  of  axial  and  transverse 
foil  strain  gages  were  used  on  front  and  back  surfaces 
in  order  to  average  out  any  through-thickness  bending. 
All  readings  were  corrected  for  gage  transverse 
sensitivity  and  strain  nonlinearity  using  the 
manufacturer's  data.  Several  specimens  were  tested 
using  two  cycles  of  tensile  loading,  unloading,  and 
reloading  at  successively  higher  maximum  strains. 

The  specific  layups  used  for  the  chracterizat ion 
phase  and  the  theoret i cal -experimental  comparisons  are 
identified  in  subsequent  sections. 

4.  LINEAR  ELASTIC  BEHAVIOR 

Elastic  constants 

The  principal  properties  are 

E , *  125.5(18.2),  v  =  0.334, 

E'=  9.38(1.36).  Gj2=  5.22(0.757)  (34) 

where  the  moduli  are  in  GPa(Msi).  Chebyschev 
polynomials  were  fit  to  all  stress-strain  curves  and  the 
first-order  coefficients  provided  the  results  in  Eq. 
(34).  Unidirectional  0‘  and  90  layups  provided 

(E,,  v)  and  E2,  respectively,  while  the  angle-ply 
laminate  (  ± 4 5 !  3^  provided  G^  from 


When  the  elastic  behavior  is  linear  Eq.  (31)  reduces  to 


G12  ‘  V2<l+Nx> 


(35) 


w  *  (Qjpj  +  G22c2  +  2','G22clt2  +  G12Y12^2  (33) 

An  alternative  characterization  in  terms  of 
stresses  using  the  dual  energy  was  also  developed.  For 
linear  elastic  behavior  it  is  equivalent  to  the  strain 
formulation.  The  elastic  nonlinearity  was  not 
Introduced  by  using  Eq.  (31)  In  Eq.  (2);  instead,  for 
simplicity  w  was  constructed  to  be  similar  In  form  to 
Eq.  (31),  btft  with  stress-dependent  functions  as  found 
from  the  uniaxial  stress  tests.  All  predictions  made 
with  this  stress  formulation  were  virtually  the  same  as 
found  from  the  strain  formulation,  Eq.  (31).  As  the 
strain  formulation  is  far  more  convenient  for  predicting 
response  of  multi-directional  fiber  laminates,  we  use 
only  this  formulation  here. 

3.  EXPERIMENTAL  PROGRAM 

The  composite  material  used  in  all  experimental 
work  was  a  graphite  fiber-reinforced  epoxy,  Hercules* 
AS4/3502  with  a  64. 2%  fiber  volume  fraction.  It  was 


where  E  and  N  are  the  Young’s  modulus  and  Poisson's 
ratio  of  the  laminate;  although  this  equation  is  the 
same  as  for  an  isotropic  material,  it  really  comes  from 
lamination  theory. 

Predictions  for  unidirectional  and  angle-ply  laminates 
Predictions  based  on  Eq.  ( 34)  and  standard  I  1  near 
theory  are  shown  by  the  continuous  lines  in  Figs.  4  and 
5. 
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Fig.  4  Modulus  and  Poisson's  ratio  of  the  angle-ply 
laminate  versus  fiber  angle 


Fig.  5  Modulus  and  Poisson's  ratio  of  the  uni¬ 
directional  laminate  versus  fiber  angle 
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Agreement  between  theory  and  experiment  Is  excellent 
except  for  the  off-axis  unidirectional  Poisson's  ratio. 

Sensitivity  study 

The  very  low  failure  strain  of  the  90"  specimen 
precludes  its  use  for  determining  Ep  as  a  function  of 
S.  In  order  to  identify  good  layups  for  finding  Ep  as 
well  as  G^p,  the  sensitivity  of  Young's  modulus  Ex  and 
Poisson's  ratio  N  to  £2  and  Gjo  was  calculated. 
Figures  6  and  7  snow  the  results  For  angle-ply 


Fig.  6  Effect  of  changes  in  G , 2  and  ^2  on  Youn9's 
modulus  Ex  and  Poisson^s  ratio  Nx  for  the 
angle-ply  laminate 


Fig.  7  Effect  of  changes  In  Gip  and  E2  on  Young's 
modulus  £  and  Poi$son*s  ratio  Nx  for  the 
unidirectional  laminate 


(u)  and  unidirectional  («)  laminates,  where  6  is  the 
angle  between  the  loading  and  fiber  directions.  1  he 
sensitivitites  are  logarithmic  derivatives;  for  example 
£x : E 2  is  the  plot  of 

a  log  E  E?  aE 

-x - r1  =  r1  TT  ( 36) 

alog  E2  Ex  3 E 2 

Thus,  an  ordinate  value  of  0.75  implies  a  1%  change  in 
Ep  produces  a  0.75%  change  in  Ex.  The  modulus  and 
Poisson's  ratio  of  the  (±45")  layup  are  practically 
Independent  of  Ep,  and  therefore  some  other  layup  is 
needed  to  obtain  Ep.  Although  the  sensitivity  to  Ep  is 
good  at  large  angles,  the  failure  strain  is  quite  low. 
The  (±30')  layup  provides  adequate  sensitivity  and  a 
relatively  high  failure  strain,  and  therefore  it  was 
selected  to  obtain  Ep.  Sensitivities  in  Fig.  7  are 
quite  good  for  most  unidirectional  layups.  The  high 
failure  strain  we  found  for  the  15"  layup  seems  to  make 
this  a  good  choice;  however,  as  discussed  later,  there 
appears  to  be  some  difference  between  the  intrinsic 
nonlinear  behavior  of  unidirectional  and  angle-ply 
laminates. 

5.  NONLINEAR  BEHAVIOR 

Haterial  characterlzat ion 

Figures  8  and  9  show  the  experimental  mechanical 


Fig.  8  Behavior  of  [Oljp  laminate 


Fig.  5  Behavior  of  [90|12  laminate 


6 


behavior  (out  to  the  ultimate  strain)  for  the  0°  anu  3C° 
specimens,  respectively.  Both  modulus  and  Poisson's 
ratio  are  seen  to  vary  somewhat  with  strain.  Data  from 
load-unload-reload  tests  showed  negligible  departures 
from  the  single-load  CT  curves,  as  noted  previously. 
Roughly  75%  of  the  nonlinear  behavior  in  the  90°  curve 
was  retraced,  and  thus  we  used  for  £gQ(c?)  In  Eq.  (29)  a 
secant  modulus  with  25%  less  nonlinearity  than  found 
directly  from  the  90"  test  data.  Second  order 
Chebyschev  polynomials  were  found  to  represent  the  0° 
and  90°  data  very  accurately;  the  uncorrected  data  for 
Poisson's  ratio  in  Fig.  9  were  adjusted  by  a  constant 
factor  so  that  the  initial  value  of  v2,  agreed  with  the 
first  expression  in  Eq.  (26c);  c  apparently,  the 
manufacturer's  transverse  sensitivity  correction  factor 
for  the  transverse  strain  gage  is  not  sufficient  to 
fully  account  for  this  sensitivity  when  the  Poisson's 
ratio  is  very  small.  Tne  polynomials  were  used  to 
evaluate  the  secant  moduli  and  tangent  Poisson's  ratio 
in  Eq.  (29).  The  top  dashed  line  in  Fig.  9  is  the 
predicted  Poisson's  ratio  from  Eq.  (27)  when  all 
nonlinearity  is  elastic. 

In  order  to  complete  the  characterization  needed 
for  Eq.  (27),  E2  =  Eo(S)  and  Gj2  =  GijlS)  were  obtained 
from  the  (t30")  and  (+45  )  data  in  Figs.  10  and  11. 


Fig.  10  Behavior  of  1+30)35  laminate 


Fig.  11  Behavior  of  1 ±4 5 ] 35  laminate 


In  Eq.  (27),  Qjj  -  Ej(o)/D  and  Q22  -  E2(S)/D,  where,  as 
an  excellent  approximation,  we  nave  used  the  linear 
elastic  constants  from  Eq.  (34)  to  find  0  =  0.992. 


After  completing  tne  cnaracterizat ion,  we  predicted 
the  inelastic  behavior  (S<0)  of  the  (  +  45)  specimen  using 
the  linear'  elastic  (fj  =  f2  -  fjo  =  1)  and  nonlinear 

elastic  representations.  No  significant  difference  was 
found.  Consequently,  one  may  use  Eq.  (35)  to  obtain  Gj2 
from  the  datainFig.il. 

To  find  E2,  Eq.  (27)  was  rewritten  in  terms  of 
specimen  coordinates  (x,y)  using  second  order  tensor 
transformations  for  the  stresses  and  strains.  Then,  py 
using  the  condition  that  y  «  0  for  an  angle-ply 

laminate,  the  axial  modulus  E^  and  Poisson's  ratio 
N  were  expressed  as  functions  of  G^2  and  E2.  This 

relult  enabled  us  to  find  G,2  and  E2  as  a  function  of 
the  axial  strain.  Although  the  (+30")  specimen  provides 
both  Gj2  and  E2,  we  found  that  the  higher  ultimate 

strain  of  the  ( ±  45 “ )  specimen  enabled  Gj2  to  be 

determined  out  to  a  17%  larger  value  of  S,  and  thus  only 
the  latter  Gi2  was  used  subsequently.  The  modulus  E2 
was  predicted  out  to  the  same  maximum  S  using  the 

polynomial . 

The  procedure  discussed  thus  far  provides  E2  and 

G 1 2  as  a  function  of  axial  strain  from  each  of  the 

specimens.  It  remains  to  relate  these  strains  to  the 

structure-change  work  S,  as  given  by  Eq.  (24).  When  the 
elastic  nonlinearity  is  neglected,  then  w  =  0  c  /2,  so 
that  this  work  for  unidirectional  and  Sn§ 1 e -p  1  y 

laminates  at  any  given  axial  strain  c  is  the  shaded 

area  Illustrated  in  Fig.  3.  With  elastic  nonlinearity 
then  Eq.  (31),  expressed  In  specimen  coordinates, 
provides  w.  In  the  neighborhood  of  the  yiitial  state 
(c  *  S  =  0),  both  wj  and  -  w  vary  as  c  ,  and  their 

difference,  S,  varies  as  Moreover,  bofti  G^  and  f2 

have  non-zero  first  order  coefficients  in  cf. 

Therefore,  in  order  to  fit  Gj2  and  E2  with  Chet.,  . 
polynomials,  S„  ;  S1  ,  rather  than  S,  was  used  as  the 

expansion  parameter. 

The  continuous  lines  in  Fig.  12  for  Gj2  and  E2  were 


Fig.  12  Transverse  and  shear  moduli  versus  cube  root  of 
the  structure-change  work 


plotted  using  a  sixth-order  polynomial  fit  in  Sp  to  the 
data  from  (+45°)  and  (+30°)  specimens;  the  units  used 
for  S  are  Msi.  The  data  points  show  moduli  found  from 
15",  30",  and  45°  unidirectional  specimens,  neglecting 
the  elastic  nonlinearity.  The  effect  of  this 
nonlinearity  was  estimated  as  being  too  small  to  account 
for  the  differences  shown  for  E2. 

Predictions 

All  material  functions  (fj,  fp,  fj2)  and  (Gj2,  E2) 
have  been  expressed  as  polynomials  in  strain  and  SR| 
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respectively.  In  order  to  predict  the  Instantaneous 
value  of  S  in  each  ply  of  a  laminate  we  may  use  (23). 
This  expression  is  a  nonlinear  algebraic  equation  for  S, 
given  the  strains.  The  strain  energy  function  is  in  £q. 
(31)  or  (33).  Prediction  of  ply  strains,  for  each  S, 
was  done  using  lamination  theory  (8).  The  numerical 
method  employed  consisted  of  solving  simultaneously  for 
the  unknown  ply  strains  and  S,  given  the  axial  strain, 
using  the  Newton-Raphson  method;  usually  only  a  few 
Iterations  (2  to  5)  .were  needed  for  convergence  to  a 
relative  error  of  10  .  We  should  add  that  since  SR  was 
used  as  the  polynomial  expansion  parameter,  we  changed 
Eq.  (23)  to  the  form 


1 

S 


R 


/  3w 

(’s« 


=  o 


(37) 


and  then  used  the  Newton-Raphson  method  to  drive  the 
left  side  to  zero;  the  indicated  division  by  SR 
Increased  the  rate  of  convergence. 

All  predictions  that  follow  are  based  on  the  Gjp 
and  Ep  in  Fig.  12  found  from  the  angle-ply  specimens; 
the  nonlinear  elastic  behavior  was  accounted  for.  When 
both  Gip  and  Ep  were  found  from  the  (iTO")  data  in  Fig. 
10,  differences  between  theory  and  experiment  in  Fig.  10 
could  not  be  discerned  whether  nonlinear  elasticity 
("nonlinear  fibers")  or  linear  elasticity  ("linear 
fibers")  was  used,  as  expected.  However,  for  the 
results  reported  here  only  one  function  was  used  from 
each  of  the  two  laminates.  Thus,  the  predictions  in 
Figs.  10  and  11  provide  a  partial  check  on  the  theory. 

Figures  13-19  show  additional  predictions 
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Fig.  13  Behavior  of  (1I5I35  laminate 


Fig.  15  Behavior  of  115]^  laminate 


Fig.  16  Behavior  of  130!^  laminate 
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Fig.  14 


Behavior  of 


( i60 1 3S  laminate 


Fig.  17  Behavior  of  MSIjp  laminate 


Fig.  18  Behavior  of  |  -t45/-45/+45/45/0/+45/-45  1 3 
laminate. 
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fig.  19  Behavior  of  | +45/-45/+45/-45/90/+45/-45 |$ 
laminate 


for  angle-ply,  unidirectional,  and  tridirectional 
laminates.  In  most  cases  there  was  little  difference 
between  the  nonlinear  fiber  and  linear  fiber  (fj=  f2= 

-  1)  predictions,  and  so  only  the  nonlinear  case  is 
shown.  Inclusion  of  the  elastic  nonlinearity  always 
improved  the  prediction.  The  "linear  variation"  in  Fig. 
18  uses  the  initial  modulus;  it  emphasizes  the  softening 
effect  of  S  in  the  45"  plies.  The  axial  stress 
prediction  using  nonlinear  fibers  is  essentially  the 
same  as  the  measured  stress. 

figure  19  shows  results  from  three  replicate  tests 
of  a  laminate  that  was  selected  for  its  susceptibility 
to  edge  delamination.  Indeed,  delamination  was  observed 
to  gradually  develop  at  the  edge  and  then  grow  inward, 
beginning  at  a  strain  of  atout  0.008.  Full  delamination 
and  specimen  failure  occurred  at  approximately  twice 
this  strain.  Transverse  cracking  (TC)  in  the  90"  plies 
initiated  just  before  edge  delaminations  were 
observed.  The  Poisson's  ratio  was  not  affected  by  the 
delamination  until  it  reached  the  centrally  located 
strain  gages. 

For  the  laminate  in  Fig.  18  the  0"  plies  are  under 
a  compressive  transverse  stress,  while  the  90"  piios  for 
that  in  Fig.  19  have  a  compressive  stress  parallel  to 
their  fibers.  In  both  cases,  this  compression  is  due  to 
the  relatively  high  Poisson's  ratio  of  the  (t45)  ply- 
pair  (cf.  Fig.  11).  No  adjustment  was  made  to  the 
nonlinear  elastic  coefficients  fj,  f2,  and  fj2  To 
account  for  this  compression,  even  though  they  were 
determined  under  tensile  loading;  as  these  coefficients 
turned  out  to  be  linear  in  strain,  a  stiffening 
nonlinearity  in  tension  becomes  a  softening  nonlinearity 
in  compression  and  vice-versa.  When  elastic 

nonlinearity  was  neglected,  the  agreement  with 
experimental  data  was  not  quite  as  good  as  when  It  was 
used. 

Delawination  analysis 

Prediction  of  the  delamination  in  the  laminate  of 
Fig.  19  may  be  accomplished  using  the  same  method  as  for 
an  elastic  material,  except  the  work  potential  density 
wj  replaces  strain  energy  density  w  (4,6).  The  two 
lowest  dotted  curves  in  Fig.  19  are  strain  energy 
release  rates  x  10  (referred  to  the  left  axis)  for  an 
edge  delamination  in  the  outer  interface  between  the  -45 
and  90  plies.  They  were  found  by  adapting  0'8r1en's 

method  (9)  to  the  present  formulation,  as  noted  above. 
Basically,  one  first  adds  the  work  potential  density 
wj1'  of  the  central,  undelaminated  section  to,-those  of 
the  two  separated  laminates,  wj/'1  and  wt  •  after 
multiplying  them  by  their  respective  volumes.  The 
delamination  work  is  then  added  to  obtain  the  total  work 


potential , 


WT=  2{wj^  (B-aJh^U  |wj2M2)  +  wj3V3ha|L  *  2GcaL(38) 


where  B  is  the  specimen  half-width,  a  is  the 
delamination  depth,  L  is  the  (assumed)  constant 
delamination  length,  h^1^  are  the  section  thicknesses 
and  Gc  is  the  delamination  fracture  energy/area.  The 
growth  condition  is  aw^/sa  =  0,  which  may  be  written  as 


G 


(39) 


where  G  is  the  energy  release  rate, 


G  =  w^V1) 


(2).  (2)  (3). (3) 


(40) 


This  expression  was  evaluated  for  two  cases  gukled  by 
the  data:  (i)  the  90°  plies  are  Included  in  wj  '  (for 
the  inner  ply-group)  and  (ii)  the  90‘  plies  are 
omitted.  The  second  case  intersects  the  horizontal  Gc 
llne  (G  =  1  lb/in  from  double  cantilever  beam  tests)  at 
e  =■  0.007,  while  the  first  one  does  so  at  t  =  0.013, 
al  shown  in  Fig.  19.  Due  to  transverse  cra/ks  and  the 
axial  fiber  compression,  the  90°  plies  may  lose  their 
ability  to  resist  transverse  contraction  (possibly  due 
to  fiber  buckling)  thus  resulting  in  early 

delamination.  Figure  19  shows  that  most  of  the  observed 
delamination  takes  place  between  these  predicted  values. 


Thermal  strains  and  Fig.  12 

Thermal  strains  due  to  cool-down  from  the  cure  at 
350F  have  been  neglected.  We  may  introduce  them  through 
the  b;,  Eq.  (16),  which  can  be  interpreted  as  free 
thermal  expansion  strains.  Clearly,  by  using  c  -  b.  in 
place  of  £  ■  in  Eqs.  (25)  -  (33),  the  effect  of'thefmal 
strains  may  be  taken  into  account.  It  should  be  noted 
that  the  omission  of  them  probably  does  not  explain  the 
differences  in  E2  in  Fig.  12;  these  differences  may 
really  be  due  to  the  effect  of  the  constraint  from 
adjacent  plies  in  reducing  microdamage  in  the  angle-ply 
laminate.  It  is  expected  that  residual  thermal  strains 
in  the  angle-ply  laminate  would  lead  to  a  lower  Ep, 
rather  than  a  higher  value,  compared  to  that  for  the 
unidirectional  laminates. 


CONCLUSIONS 

In  the  simplest  form  of  the  theory  used  here,  the 
unidirectional  ply  was  modeled  as  a  linear  elastic 
material  with  principal  shear  and  transverse  moduli  that 
vary  with  one  scalar  parameter  (which  Is  equal  to  the 
applied  work  less  the  strain  energy  for  a  unit 
volume).  A  modified  form  of  this  theory  was  introduced 
to  account  for  elastic  nonlinearity  arising  from 
Straining  parallel  and  perpendicular  to  the  fibers 
(possibly  due  to  initial  fiber  waviness).  Predicted  and 
experimental  laminate  results  for  various  layups  were 
found  to  be  in  good  agreement.  In  order  to  predict 
delamination,  it  was  observed  that  conventional  energy 
release  rate  analysis  may  be  used,  but  the  work  density 
function  replaces  strain  energy  density.  Although 
predictions  of  unloading  and  reloading  response  were  not 
made,  a  possible  approach  was  described.  Viscoelastic 
behavior  was  neglected;  relatively  simple  approaches 
have  been  proposed  elsewhere  to  account  for  this 
behavior  at  crack  tips  (4)  and  in  the  continuum  (10). 
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